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ABSTRACT 

We present high S/N echelle spectra of the recently discovered ultraluminous 
QSO APM 08279+5255 and use these data to re-examine the abundance 
of Carbon in Lya forest clouds. In agreement with previous work, we find 
that approximately 50% of Lya clouds with hydrogen column densities 
log iV(H I)> 14.5 have associated weak C IV absorption with log N(C IV) £ 12, 
and derive a median N(C IV)/iV(H I)= 1.4 x 10~ 3 . The agreement with earlier 
estimates of this ratio may be somewhat fortuitous, however, because we show 
that previous analyses have probably overestimated the number of Lya clouds 
which should be included in this statistic. 

We then investigate whether there is any C IV absorption associated with 
weaker H I column densities by stacking 51 C IV regions corresponding to 
51 Lya lines with 13.5 < logiV(H I) < 14.0. The co-added spectrum has 
S/N ~ 580 but shows no composite C IV absorption. In order to understand 
the significance of this non-detection we have stacked together 51 theoretical 
C IV A1548 lines with individual values of column density and velocity dispersion 
scaled appropriately from the values measured in the corresponding Lya lines. 
We find that even if the typical value N(C W)/N(R I)= 1.4 x 10" 3 applies to 
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these lower column density clouds, the corresponding signal in the stacked C IV 
region is smeared by the likely random difference in redshift between C IV and 
Lya absorption and becomes very difficult to recognize. This seems to be a 
fundamental limitation of the stacking method which may well explain why in 
the past it has led to underestimates of the metallicity of the Lya forest. 

We also analyze our spectra with the pixel-by-pixel optical depth technique 
recently developed by Cowie & Songaila (1998) and find evidence for net C IV 
absorption in Lya clouds with optical depths as low as r(Lya) = 0.5 — 2, as 
these authors did. However, we show with simulations that even this method 
requires higher sensitivities than reached up to now to be confident that the 
ratio iV(C IV)/iV(H I) remains constant down to column densities below 
log iV(H I)~ 14.0 . We conclude that the question of whether there is a uniform 
degree of metal enrichment in the Lya forest at all column densities has yet to 
be fully answered. Future progress in this area will probably require concerted 
efforts to push further the detection limit for C IV lines in selected bright QSOs. 

Subject headings: galaxies: formation - galaxies: intergalactic medium - 
quasars: absorption lines - quasars: individuals (APM 08279+5255) 



1. Introduction 

Nearly thirty years after its discovery by Lynds (1971), the Lya forest continues to 
provide fertile ground for observational and theoretical studies. New echelle spectrographs 
on large telescopes, most notably the High Resolution Spectrograph (HIRES, [Vogt et aLl 
1994|) on Keck I, have provided spectra of exceptional quality and detail. This has permitted 



extensive study of many important properties of the Lya forest as comprehensively 
reviewed by Rauch (1998). Complementary to this quality observational data, cosmological 
simulations which utilize gas hydrodynamics have improved our understanding of the nature 
of the forest clouds. These simulations have shown that the formation of the Lya forest 
is a natural consequence of the growth of structure in the universe through hierarchical 
clustering in the presence of a UV ionizing background (e.g. |Cen et al. 1994}; [Petitjean^ 



Mucket fc Kates 1995) ; [Hernquist et al. 1996| ; [Bi fc Davidsen 1997] ) . In this class of models, 



low column density H I clouds (log iV(H I)<14, where iV(H I) is measured in cm -2 ) are 
found preferentially in voids whilst stronger lines are associated with higher density clouds 
which arise in filamentary structures around collapsed objects. 



Important clues to the origin of Lya clouds can be gleaned from their clustering 
properties and chemical composition. Early investigations of the Lya forest found no 
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evidence of metals associated with the H I clouds, suggesting that they consisted of pristine 
material ( [Sargent et al. 1980| ). More concerted efforts however, did provide some evidence 
for metals in the forest either individually (Meyer & York 1987) or in a stacked spectrum 
(Lu 1991). The realization that metal enrichment is in fact widespread in high column 
density Lya clouds has come relatively recently with the availability of HIRES on the Keck 
I telescope. A number of studies (e.g. [Tytler et al. 19951 ; |Cowie et al. 1995| ; |5ongaila~Sz 



|Cowie 1996|) have found that C IV absorption is seen in approximately 50% of forest clouds 
with log iV(H I)>14.5 and in about 90% of clouds with log N(R I)>15. Evidently the 
intergalactic medium (IGM) has been enriched by the products of stellar nucleosynthesis 
even at redshifts as high as z ~ 3.5. Photoionization models (e.g. [Hellsten et al. 1997|) can 



reproduce the observed C IV/H I ratio at z ~ 3 for a typical Carbon abundance^ [C/H] = 
—2.5 . Rauch, Haehnelt, & Steinmetz (1997) predict an order of magnitude scatter in [C/H] 
in models where random sight-lines are cast through protogalactic clumps, a prediction 
corroborated by the data. 

Two scenarios have been suggested to explain the presence of metals in the Lya forest; 
early pre-enrichment by a widespread episode of star formation, possibly associated with 
Pop III stars, or in-situ enrichment whereby the H I cloud is enriched locally either by 
contamination from a nearby galaxy or by star formation within the cloud itself ( |Ostriker| 
|fc Cncdin 1996| ; Cnedin fc Ostrikcr 1997| ). Cosmological simulations of a two phase ISM 



recently performed by Gnedin (1998) confirm the results of Gnedin & Ostriker (1997) 
that the dominant process for transportation of heavy elements into the IGM is mergers 
of protogalaxies in dark matter halos. One may expect that for the scenario of in-situ 
enrichment, where metals are synthesized relatively close to the Lya clouds being observed, 
the metallicities would be highly inhomogeneous depending on the proximity of a given 
Lya cloud to a star-forming region and its past merger history. If, however, the bulk of 
the metals is formed by an episode of Pop III star formation at sufficiently high redshift 
(z ~ 14 in the models of Ostriker & Gnedin 1996), a more uniform [C/H] may have been 
established by z ~ 3, if the mixing mechanism is efficient. 

Discrimination between these two enrichment scenarios is probably best addressed in 
the low column density regime. Since the optical depth of Lya clouds is roughly indicative 
of the baryon overdensity, studying low column density lines will give us an insight into the 
low mass systems. In the simulations of Gnedin & Ostriker (1997), for example, a sudden 
drop in metallicity is predicted for Lya clouds with log iV(H I)<14 where the star formation 
rate is expected to be much lower than in higher density clouds. However, determining 
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the [C/H] abundance in weak Lya clouds has been observationally challenging due to 
the extreme weakness of C IV absorption when log N(H I) £ 14. In order to overcome 
the difficulty of a direct C IV detection, this problem has recently been tackled with two 
different approaches which have produced apparently conflicting results. 

Lu et al. (1998) tackled the problem by stacking almost 300 C IV regions to produce a 
composite spectrum, a method first applied to QSO spectra by Norris, Peterson & Hartwick 
(1983). Having found no C IV in their summed spectrum (which had a final S/N ratio 
of 1860) Lu et al. used Monte Carlo simulations to place an upper limit on the Carbon 
abundance in clouds with log iV(H I)< 14 of [C/H]< —3.5. A different conclusion was 
reached by Cowie & Songaila (1998) who used distributions of Lya optical depths, r(Lya), 
to build distributions of corresponding t(C IV) on a pixel-by-pixel basis. A comparison 
between the r(C IV) distribution and a reference 'blank' distribution showed a residual 
signal sufficiently strong to be consistent with [C/H]~ —2.5 for Lya clouds with log iV(H I) 
as low as 13.5. Thus, whilst Lu et al. interpreted their results as evidence against a uniform 
enrichment of the IGM at z ~ 3, Cowie & Songaila proposed that transportation/ejection 
mechanisms from early sites of star formation are much more efficient than anticipated, 
in order to pollute Lya clouds uniformly over a range of N(H I) of nearly 4 orders of 
magnitude. 

In this paper we aim to address the contradiction between these two earlier analyses by 
taking another look at the metallicity in the Lya forest using new data which are among 
the best ever obtained for this purpose; almost 9 hours of HIRES observations (described 
in §2) of the ultra-luminous QSO APM 08279+5255. In §3 we define the redshift interval 
over which we search for C IV absorption; the very high S/N of the spectrum in the C IV 
region allows us to refine earlier estimates of the typical C IV/H I ratio in Lya lines with 
log iV(H I)> 14.5 (§4). In §5 we investigate whether this high quality, independent data set 
can help resolve the current discrepancy between the results of Lu et al. (1998) and Cowie 
& Songaila (1998), thus shedding some light on the origin of the metals observed in low 
column density Lya clouds. We summarize our main results in §6. 



2. Observations and Data Reduction 

The target for these observations is APM 08279+5255, an ultra-luminous Broad 
Absorption Line (BAL) quasar discovered by Irwin et al. (1998) during a survey of Galactic 
halo carbon stars. The emission redshift z cm = 3.87 measured by Irwin et al. from the C IV 
A1549 and N V A1240 emission lines has recently been refined to z cm = 3.9110 by Downes 
et al. (1999) who detected CO emission. Given the difficulty in measuring a precise redshift 
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from the UV emission lines which are affected by the BAL phemomenon, the difference 
is probably not significant; in the following analysis we adopt z = 3.9110 as the systemic 
redshift of the QSO. Positionally coincident with an IRAS Faint Source Catalog object with 
a 60 micron flux of 0.51Jy, APM 08279+5255 has an optical R magnitude=15.2 and an 
inferred bolometric luminosity of 5 x 10 15 L Q making it the most luminous object currently 



known. Adaptive optics imaging with the CFHT ([Ledoux et al. 1998|) has resolved APM 
08279+5255 into two images separated by 0.35 arcsec in the NE-SW direction with an 
intensity ratio Ine/Isw — 1-2 ± 0.25 in the if -band. 

For our purposes, APM 08279+5255 is a nearly ideal background source for 
investigating the Lya forest and associated metal absorption lines at high S/N ratio and 
resolution. The data presented here were obtained with HIRES on the Keck I telescope 
on three runs in April and May 1998. Details of the observations are presented in Table 
1. APM 08279+5255 is unresolved in these observations. The cross disperser and echelle 
angles were used in a variety of settings to give almost complete wavelength coverage 
between 4400 and 9250 A. The emission spectrum of a Th-Ar hollow cathode lamp provided 
a wavelength reference and a continuum source internal to the spectrograph was used for a 
first order correction of the echelle blaze function. 

The data were reduced using Barlow's (1999, in preparation) customized HIRES 
reduction package (HAR). The individual, sky-subtracted spectra were mapped onto a 
linear wavelength scale with a dispersion of 0.04 A per wavelength bin and then co-added 
with a weight proportional to their S/N. Finally, this co-added spectrum was normalized by 
fitting a cubic spline function with STARLINK software to continuum regions deemed to 
be free of absorption. Even in the dense forest of Lya lines at wavelengths below that of 
Lya emission, continuum windows can be identified at the high resolution of our data. 

The final co-added spectrum has a resolution of 6 km s _1 , sampled with ~ 3.5 
wavelength bins, and S/N between 30 and 150. Over the range of interest for C IV 
absorption (see §3 below) the typical signal-to- noise ratio is S/N ~ 80. At the mean 
<2abs = 3.41 the 5a detection limit for the rest-frame equivalent width of C IV absorption line 
with the typical FWHM = 22 km s _1 (see later) is Wo (1548)= 3 mA, which corresponds to 
a column density log N(C IV) = 11.88 . Thus these data are comparable to those obtained 
by Songaila & Cowie (1996) for Q0014+813 and Q1422+231, their best observed QSOs. 

The full spectrum is described elsewhere (Ellison et al. 1999) and is available via 
anonymous ftp from ftp.ast.cam.ac.uk (pub/sara/APM0827). 
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3. Sample Definition 

APM 08279+5255 belongs to a class of objects whose spectra are known to exhibit 
broad, high velocity, intrinsic absorption lines. Before proceeding we must therefore 
define a working wavelength interval over which we can study C IV and Lya lines in our 
spectrum and be confident that we are not confusing them with absorption from ejected 
material. The upper limit of this interval (for C IV) is 7280 A (z max = 3.701), at the 
blue edge of the broad C IV absorption trough which corresponds to an ejection velocity 
t>ej — 13100 km s^ 1 relative to the systemic redshift. The lower limit for the C IV interval 
is 6365 A (z m i n = 3.109), the wavelength of Ly(3 emission. We excluded the region of the 
spectrum blueward of this wavelength in order to avoid confusion of Lya with Lj(3 lines. 
We further excluded a small region between 6860 and 6950 A which is contaminated by the 
atmospheric B band. 

In principle, some of the C IV doublets in our sample could still be related to the BAL 
phenomenon and be intrinsic to the QSO rather than intervening absorbers, since ejecta 
have been observed at velocities as large as 60,000 km s" 1 . However, we found no evidence 
of this in APM 08279+5255. For example, none of the C IV systems at z^ < z mSuX show 
any indication of partial coverage of the QSO (one of the signatures of intrinsic absorption), 
whereas we do see such an effect for systems at z^ s > z max . Most importantly, as we will 
discuss below, the density of C IV absorbers per unit redshift between z m \ n and z max is 
entirely compatible with that measured towards non-BAL QSOs — there is no excess of 
C IV absorbers within our working wavelength interval. 

4. C IV In Lya Clouds with log JV(H I) > 14.5 

Before addressing the question of CIV absorption in Lya clouds with 13. 5< log N(H I) 
< 14.0, the first stage of our analysis is to determine the typical N(C IV)/iV(H I) ratio in 
Lya clouds with log iV(H I) > 14.5. The aim is to provide an independent measure of this 
quantity, and therefore of the metallicity of Lya clouds, for comparison with the results of 
Cowie et al. (1995) and Songaila & Cowie (1996). 

These authors found that at a detection limit log N(C IV) > 12.0 approximately half 
of the Lya lines with log V(H I)> 14.5 have associated C IV absorption. Lya lines above 
this column density cut-off are saturated and therefore not on the linear part of the curve of 
growth. Therefore, an accurate iV(H I) cannot normally be measured without higher-order 
Lyman lines such as Lj(3 or Ly7 that are not saturated. Cowie et al. (1995) circumvented 
this difficulty by assuming that at log 7V(H I)> 14.5 the residual flux in the core of a Lya 
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absorption line is rj < 0.025 for the mean value of the Doppler parameter of b = 34 km s _1 
determined by Carswell et al. (1991)[j. This approach is of course critically sensitive to a 
precise determination of the zero level. 

Figure 1 shows the Lya forest in APM 08279+5525 between z ahs = 3.109 and 3.701; 
in this redshift interval there are 58 Ly« lines with r/ < 0.025. However a search for 
corresponding C IV absorption yielded only 22 C IV systems (38%) at a significance level 

Of > 5(7. 

This apparent discrepancy led us to examine how accurately the 77 < 0.025 cut 
actually selects Lya clouds with log iV(H I)> 14.5. To this end we used the profile fitting 
package VPFIT ( [Webb 1987|) to determine the column density, redshift and b parameter 



of Lya lines with < 0.025. We found that although these lines are saturated, in most 
cases there is sufficient information in the absorption profiles for VPFIT to converge to a 
reliable solution. Fitting all lines which do not have a flat core at zero residual intensity, 
showed that 22 out of 58 Lyct lines with r ^ < 0.025 do in fact have log N(R I)< 14.5. We 
have reproduced some examples in Figure 2 and Table 3. The revised sample then consists 
of 36 Lya lines with measured log N(H I)> 14.5; 20 of these (56%) have associated C IV 
systems. This fraction is in good agreement, although somewhat fortuitously, with the 
value of 58 ± 8% reported by Songaila & Cowie (1996). 

All of the C IV systems found within our redshift interval were fitted with Voigt profiles 
using VPFIT (all C IV lines were unsaturated, thus allowing accurate determinations of 
their column densities) and the model parameter fits for z, N(C IV), and b are listed in 
Table 2. Fig 3 is an atlas of the 23 C IV systems whose corresponding Lya lines lie within 
our redshift interval; 22 of these have Lya with rj < 0.025. Only the 20 C IV systems with 
a measured log N(H I)> 14.5 are considered in the statistics below. 

In order to derive the typical N(C 1V)/N(H I) ratio in Lya clouds with log iV(H I) 
> 14.5, we follow the procedure used by Cowie et al. (1995). For an H I column density 
distribution of the form 

n(N)dN oc N-^dN (1) 



with /3 ~ 0.7 (Petit. jean et al. 1993|) , the median JV(H I) is 2 P times the minimum value 



or, in our case, 5 x 10 14 cm 2 . Our 5er detection limit for C IVA1548 corresponds to 
N(C IV)= 7 x 10 11 cm" 2 for the median b(C IV) = 13.4 km s" 1 (see Table 2). Since at this 
detection level approximately half of the Lya clouds have associated C IV absorption, we 



7 As usual b — \f2o, where a is the one-dimensional velocity dispersion of the absorbers projected along 
the line of sight 
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deduce a median N(C IV)/iV(H I)~ 1.4 x 10 3 . For comparison, Songaila & Cowie (1996) 
reported median values of N(C IV)/iV(H I) between 1.6 x 1CT 3 and 2.8 x 10~ 3 . 

In Figure 4 we show the column density distribution of C IV lines again assumed to be 
a power law of the form 

f(N)dN = BN a dN (2) 

where f(N) is the number of systems per column density interval per unit redshift path. 
The redshift path (used instead of z in order to account for comoving distances) is given by 
X(z) = + z ) 2 ~ 1] f° r Qo — 0.0 A maximum likelihood (e.g. Schechter & Press 1976) 
fit to the total (unbinned) sample of 20 C IV systems yielded a = —1.0 ± 0.1 for column 
densities in the range 11.74 <log N(C IV) < 14.25. This value of a, which is shown as 
a solid line in Figure 4, would indicate a flatter power law distribution than the value 
a = —1.5 (dashed line in Figure 4) deduced by Songaila (1997) from her analysis of a larger 
sample of 81 C IV absorbers towards seven QSOs. Formally, the difference is significant. A 
Kolmogorov-Smirnov (K-S) test returns a 76% probability that our data points are drawn 
from a distribution with a = —1.0, but only a < 1% probability that they arise by chance 
from a distribution with a = —1.5. However, there are significant variations in the statistics 
of C IV absorbers between different sight-lines in Songaila's sample and the integral of our 
power-law distribution, J2X(C IV) /AX = 2.3 x 10 14 cm -2 , is the same as that reported 
by Songaila for one of the sight-lines in her study, towards Q0956+122. Thus the difference 
between our best fitting value of a and that determined by Songaila (1997) may just be due 
to the limited statistics of the current samples. 



5. Lya Clouds with 13.5 < log JV(H I)< 14.0 

5.1. The stacking method 

We now turn to the question of whether Lya clouds with log iV(n I)< 14.0 have 
a significantly lower metallicity than those with log iV(n I)> 14.5. Lu et al. (1998) 
tackled this problem by stacking together the C IV regions corresponding to ~ 300 Lya 
lines towards nine QSOs; no significant signal was found in the final composite spectrum. 
In constructing their sample, Lu et al. targeted column densities in the range 13.5 < 
log N(H I) < 14.0 and, for expediency, assumed that this range corresponds to values of 
residual flux 0.05 < r/ < 0.38. Although our statistics are more limited, we can improve on 
previous analyses by constructing our sample more carefully and by simulating the results 



qo=0 and Ho = 65 km s 1 Mpc 1 are adopted throughout this paper unless otherwise stated. 
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of the stacking process for our observed distributions of column densities and 6-values, as 
we now discuss. 

We started by using VPFIT to fit all unsaturated Lya lines in the spectrum of 
APM 08279+5525 between z ahs = 3.109 and 3.701 . This produced a sample of 86 lines 
with 13.5 <log N(H I) < 14.0. Had we selected on the basis of rj rather than iV(H I), 
approximately 30% of lines in the desired column density range would have been missed, 
and 35% of lines outside the range would have been erroneously included. The lines missed 
are mainly lines with large values of b, and lines which are partially blended but for which 
VPFIT can still satisfactorily recover the individual values of 7V(H I). In the next step we 
visually inspected the C IV region corresponding to each Lya line to assess whether it is 
suitable for stacking; regions contaminated by atmospheric features or blended with other 
lines were excluded. This produced a final list of 51 C IVA1548 regions and 40 C IV A 1550 
regions each 7 A wide, which could be stacked. After reducing each portion of the spectrum 
to the rest frame using the values of z a t, s appropriate to each Lja line (as measured with 
VPFIT), the regions were summed and renormalized. The final co-added spectrum has S/N 
= 580 and is reproduced in Figure 5a. 

It is clear from Figure 5a that no absorption is detected in the stacked spectrum. In 
order to assess the significance of this non-detection, we compare the stacked spectrum with 
that obtained by adding together the same number of C IV lines with strengths predicted 
for different values of N(C IV)/7V(H I). In performing these simulations we again made use 
of the available information on the distribution of column densities and 6-values appropriate 
to each line, rather than assuming representative values of these quantities, as was done in 
earlier analyses. Specifically, we simulated the absorption profiles of 51 C IVA1548 lines 
assigning to each a column density N(C IV)j = N(H I)j x A, where 7V(H I); is the neutral 
hydrogen column density returned by VPFIT for the ith Lya line, and A is the adopted 
ratio N(C IV)/7V(H I). We repeated the simulations for three values of A = 2.5 x 10~ 3 
(Cowie & Songaila 1998), 1.4 x 10~ 3 , the median value deduced here for clouds with 
log iV(H I)> 14.5, and 0.7 x 10~ 3 , a factor of two lower than our median. 

Concerning the value of b to be assigned to each C IV line in the simulations, we could 
in principle consider two possibilities. If b reflected primarily large-scale motions of the 
absorbing gas, b(C IV) i = b(H /)«. On the other hand, if b has mainly a thermal origin, the 
6-values would scale as the square root of the atomic mass and b(C IV) i = \ j\J~Y2 x b{H /)«. 
The real situation will be somewhere between these two extremes. Rauch et al. (1996) 
found mean and median values of b(C IV) near 10 km s _1 from their analysis of 208 
C IV absorption components, and concluded that both thermal motions and bulk motions 
contribute to the line broadening. Our distribution of b(C IV) for Lya clouds with 
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log iV(FI I)> 14.5 has a median value b(C IV) = 13.4 km s -1 whereas for 13.5 < 
log N(R I)< 14.0 we find a median b(H I) = 28.7. In the simulations we therefore adopted 
b(C IV") i ~ 1/2 b(H I)i, with the implicit assumption that there is no significant change in 
the distribution of 6-values with column density. 

Theoretical C IVA1548 absorption line profiles were computed in this manner for the 
51 Lya lines in our sample, convolved with the instrumental resolution, and then co-added. 
The resulting composite profiles, degraded with random noise corresponding to the S/N = 
580 of the real stacked data, are shown in Figures 5b, c, and d. Also shown in these panels 
is the corresponding theoretical absorption profile for the mean N(C IV). 

The case A = 2.5 x 10~ 3 (Figure 5b) produces a composite absorption feature which is 
significant at the 7a level; such a high value of the N(C IV)/iV(H I) ratio would therefore 
appear to be excluded by our observations. At these low optical depths the equivalent width 
of the co-added C IVA1548 line scales linearly with A. Thus, if the median A = 1.4 x 10 -3 
found above for Lya clouds with log N(R I) > 14.5 also applied at log N(R I) = 13.5 — 14.0, 
we may also expect a detectable signal (4<j, Figure 5c), whereas no composite absorption 
would be recognized if A = 0.7 x 10~ 3 (Figure 5d). 

In the above simulations we have assumed that z(C IV) i = z(H J)j. However, in 
reality there is likely to be a dispersion of values of Azi = z(C IV) i — z(H J)j. We have 
already seen that b(C IV) i ~ 1/2 b(H I) i: indicative of the fact that Lya absorption takes 
place over a wider velocity range than C IV. If we compare the measured values of z(C IV) i 
and z(H I)i in the 20 Lya clouds with log N(R I) > 14.5 which show C IV absorption (§4 
above), we find a distribution of values of Azi with a z = 4 x 10~ 4 which corresponds to a 
velocity dispersion of 27 km s" 1 . We therefore repeated the simulations above assuming 
that the C IVA1548 lines to be stacked are drawn at random from a Gaussian distribution 
of Azi with this value of a; the results are reproduced in Figures 6b, c, and d. 

Clearly, once this redshift dispersion is introduced in the sample, the signal in the 
stacked spectrum is blurred to the point where it is questionable whether it would be 
detected even if N(C TV)/N(H. I) = 2.5 x 10~ 3 . Although the average residual intensity in 
Figure 6b is less than one, the absorption is so diffuse that it becomes difficult to distinguish 
it from small fluctuations of the continuum level. We thus seem to have identified a 
fundamental limitation of the stacking method for the detection of very weak absorption 
features. Unless it can be shown that Az decreases with log iV(H I), it would appear 
that the uncertainty in the exact redshift of the metal lines relative to Lya casts serious 
doubt on the interpretation of non-detections in composite spectra, even at a signal-to-noise 
ratio well in excess of that of the observations reported here. In any case, this probably 
explains why the first attempts at stacking QSO spectra to search for weak C IV lines (e.g. 



- 11 - 



Lu 1991) underestimated the typical iV(C YV)/N(H I) ratio subsequently measured with 
HIRES spectroscopy. It is interesting to note in this context that the only case where the 
stacking method has produced an incontrovertible detection is when Barlow & Tytler (1998) 
co-added HST FOS spectra to detect C IV associated with Lja clouds at low redshift 
((z) ~ 0.5). At the coarse resolution of the FOS spectra — a factor of > 25 worse than that 
of HIRES data — the velocity difference between C IV and Lja becomes a secondary effect. 



5.2. Pixel-by-pixel Comparison 

Cowie & Songaila (1998) have recently proposed a novel way to look for weak 
absorption features associated with Lya forest lines. The method involves constructing 
cumulative distributions of optical depths for all the pixels in a spectrum where (in this 
case) C IV A1548 absorption may be found, that is all pixels at wavelengths 

1548.195 , 

j = 121^670 X A * (3) 

where Aj is the wavelength of the ith pixel in the Lja forest. Since no C IV absorption is 
expected for Lja pixels with optical depth r(Lja) < 0.1, the distribution of corresponding 
values of r(C IV) provides a reference 'blank' sample. One can then examine the difference 
between the distribution of r(C IV) for a particular range of r(Lya) values and the blank 
sample to determine whether a residual signal is present and indeed Cowie and Songaila 
(1998) reported excess C IV absorption for all Lja optical depths, from r(Lya;) = 20 to 0.5 . 

We followed closely the procedure outlined by Cowie & Songaila (1998) in applying 
their method the spectrum of APM 08279+5525. The results are reproduced in Figure 7, 
where the left-hand panel shows the distributions of r(C IV) for the same three ranges of 
r(Lyo;) considered by Cowie & Songaila. The differential distributions relative to the blank 
sample are shown in the right-hand panel; here the deficit of pixels with negative values of 
t(C IV) and the excess of pixels with positive r(C IV) give an 'S wave' pattern which is 
indicative of C IV absorption associated with Lja forest lines. As can be seen from the 
Figure, we confirm the finding by Cowie & Songaila of a C IV signal in all three intervals of 
r(Lyo;) considered. 

It is of interest, of course, to ask how changes in the A^(C IV)/iV(H I) ratio would be 
reflected in the optical depth distributions shown in Figure 7. We are in a good position to 
explore this question having quantified the absorption parameters of the entire Lja forest 
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in APM 08279+55250. To this end, we simulated the C IV absorption associated with 375 
Lya lines for which VPFIT returned values of iV(H I), z, and b, and applied Cowie and 
Songaila's technique to determine the distributions of r(Lya) and r(C IV). Since the aim 
is to test whether there is a sudden change in the metallicity of Lya clouds at low optical 
depths, we considered two possibilities: 

(1) N(C TV)/N(H I)= 2.5 x 1(T 3 for all values of iV(H I); and 

(2) a drop in iV(C IV)/iV(H I) by a factor of 10 for JV(H I)< 14.0 . 

We carried out the simulations for two cases which we term 'ideal' and 'realistic'. In the 
ideal case, we did not include the effects of noise (that is, we assumed infinitely high 
S/N) and we adopted b(C IV)i = b(H J)j and z(C IV) i = z(H J)j . In the realistic 
case, we introduced random noise in the simulations to reproduce the S/N = 80 typical 
of our HIRES data. We also assumed b(C IV) i = 1/2 b(H I)i and a dispersion of values 
Azi = z{C IV)i - z(H I)i with a z = 4 x 10~ 4 as found for Lya lines with log N(E I) > 14.5 
(§5.1) . The results of these simulations are reproduced in Figure 8. Since the assumed drop 
in the N(C TV)/N(H I) ratio at N(R I)< 14.0 (case 2) would affect primarily pixels with 
-r(Lya) £ 1, for clarity we show the results at low optical depths in two subsets: r(Lya) 
< 1 and 1 < r(Lya) < 2 . 

The top panel in Figure 8 shows that in the 'ideal' case the optical depth method 
developed by Cowie and Songaila (1998) would indeed be sensitive to changes in the Carbon 
abundance of Lya clouds (open squares). Under more realistic conditions (bottom panel), 
however, the break below r(Lya)~ 1 is softened, primarily by the effect of noise which 
moves pixels with 1 < r(Lya) < 2 into the lower optical depth interval and vice versa. 
Furthermore, at the typical S/N of the present data, the error associated with the lowest 
t(C IV) bin is so large that the two cases considered (filled dot and open square) can no 
longer be distinguished with confidence. Although the data (open star) apparently favour a 
constant N(C IV)/iV(H I), we feel that the existence of a break below r(Lya)~ 1 can only 
be tested reliably with higher S/N observations. 

6. Summary and Conclusions 

We have presented a high resolution (~ 6 km s _1 ) and high S/N (~ 80) spectrum 
of the ultraluminous BAL QSO APM 08279+5255 obtained with HIRES on the Keck I 
telescope. These data, which have a sensitivity to C IV AA1548, 1550 lines with rest frame 



'Within the redshift limits designed to exclude BAL material, as discussed at §3 above 
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equivalent widths as low as ~3 mA (5a), have been analyzed with a view to reassessing the 
metallicity of the Lya forest at z ~ 3.4. Our principal results are as follows. 

In agreement with previous analyses, we find that approximately 50% of Lya clouds 
with log iV(H I)> 14.5 have associated C IV systems with log N(C IV) £ 12.0; we deduce a 
median N(C IV)/AT(H I)~ 1.4 x 10~ 3 . However, we point out that the criterion previously 
adopted to identify such Lya clouds — a residual flux in the line core rj < 0.025 — is an 
inadequate approximation which can miss a significant proportion of the lines. Even 
though the absorptions are close to saturation, profile fitting methods which simultaneously 
determine the column density and velocity dispersion of the absorbers are preferable for 
a clean definition of the sample (for data of sufficiently high S/N and resolution, such as 
those presented here). 

We have stacked the C IVA1548 regions corresponding to 51 Lya lines with 13.5 < 
log N(H I) < 14.0 but find no detectable signal in the composite spectrum. In order to 
understand the significance of this null result, we have performed simulations in which 
we stack 51 synthetic C IVA1548 lines with absorption parameters resembling as closely 
as possible those of the lines we are trying to detect. Specifically, each C IV line to be 
stacked was assigned values of column density and velocity dispersion scaled directly from 
those of the corresponding Lya line. These simulations show that we should be able 
to detect a marginally significant (4a) absorption feature in the co-added spectrum if 
N(C IV)/iV(H I)= 1.4 x 10" 3 , as in clouds with log iV(H I) > 14.5, provided that there is 
no random difference Az between the redshifts of C IV and Lya. 

However, when such a difference is introduced in the simulations the co-added signal is 
effectively washed out to the point where it can no longer be recognized. This is the case 
even for a relatively small dispersion of values of Az (a = 27 km s _1 ) which is entirely 
plausible given the observed distribution of Az in the higher column density clouds where 
C IVis detected. We propose that this uncertainty in the exact registration of weak signals 
is a serious problem which ultimately limits the usefulness of stacking procedures. 

We also analysed our data with the pixel-by-pixel optical depth method recently 
developed by Cowie & Songaila (1998) and, in agreement with these authors, we do 
find a signal indicative of C IV absorption in Lya clouds with optical depths as low as 
r(Lya) = 0.5 — 2 . However, our simulations show that a higher S/N than achieved up to 
now is required to detect with confidence even a marked break in the N(C TV)/N(R I) 
ratio below log JV(H I) < 14.0 . 

We conclude that the question of whether the abundance of Carbon in the Lya forest 
is uniform at all column densities has yet to be fully answered. Probably this question is 
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best addressed by pushing the detection limit for C IV absorption lines below the limits 
reached so far. Thanks to the extraordinary luminosity of APM 0827+5255 this approach 
is now feasible with a concerted effort of HIRES observations. 
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Fig. 1 — The Lya forest in APM 08279+5255 within the redshift range z = 3.109 - 3.701 
of our analysis. The y-axis is normalized counts. Vertical tick marks near the continuum 
level indicate Lya lines with residual flux rj < 0.025; solid and broken tick marks are used 
to indicate lines with and without associated C IV AA1548, 1550 absorption, respectively. 
Solid triangles flag Lya lines with column density log iV(H I)< 14.5 even though r/ < 0.025 
(see text). The horizontal black line between 5387 and 5456 A indicates the region of the 
forest omitted from the analysis because the corresponding C IV wavelengths fall within the 
atmospheric B band. 
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Fig. 1. — Continued 
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Fig. 2. — Examples of Lya absorption lines with rf < 0.025 and log iV(H I)< 14.5. 
Theoretical line profiles computed with VPFIT are shown as a continuous line in each panel. 
Parameters of the fits are collected in Table 2. 
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Fig. 3. — Atlas of C IV doublets associated with Lya lines in our working region. The x-axis 
is wavelength in A and the y-axis is normalized counts. Thin dotted lines show the profile 
fits with the parameters listed in Table 2. The weakest C IV systems are indicated with tick 
marks to guide the eye. 
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Fig. 3. — Continued 
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Fig. 3. — Continued 
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Fig. 4.— Column density distribution of C IV absorbers in APM 08279+5255. The data 
have been grouped into bins of width 0.5 in log N(C IV) for display purposes only. The 
top panel shows the cumulative distribution and the bottom panel is the differential. A 
maximum likelihood analysis of the distribution, assumed to be a power-law of the form 
f(N)dN = BN a dN, returned a best fit index a = —1.0 indicated by the solid line. The 
dashed line in the bottom panel (a = —1.5) is the slope reported by Songaila (1997) for her 
larger sample. 
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Fig. 5. — (a): Co-added spectrum obtained by stacking 51 C IV regions corresponding to Lja 
lines with 13.5 < log JV(H I) < 14.0. The resulting S/N is 580. (b), (c), and (d): Simulated 
stacked spectra for different values of N(C IV)/AT(H I), as indicated. The continuous line 
in each panel shows the noise-free absorption profile for a single C IV line with the mean 
N(C IV). These simulations assumed no redshift difference between C IV and Lya lines. 



-26- 



Stacked Spectra 



1.01 



1-00 - „ - 



0.99 



0.98 




Simulated Spectra 



1.01 



1.00 



0.99 



0.98 



C IV/H I = 2.5 xlO" 3 



(b) 
i 



1548 



1548.5 



1548 



1548.5 



Simulated Spectra 



1.01 



1.00 MiLTu 



0.99 - 



0.98 



C IV/H I = 1.4 xlO" 3 



(c) 
l 



1548 1548.5 
Wavelength (A) 



Simulated Spectra 



1.01 
1.00 
0.99 
0.98 




C IV/H I = 0.7 xlO" 3 



(d) 



1548 1548.5 
Wavelength (A) 



Fig. 6. — As Figure 5, except that the simulated spectra include allowance for a random 
redshift difference Az between the centroids of C IV and Lya absorption. Values of Az were 
drawn from a Gaussian distribution with a = 27 km s _1 . 
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Fig. 7. — Left: Distribution of C IV optical depths for different ranges of Lya optical depth, 
as indicated. Right: Differential distributions of C IV optical depths, relative to that for 
-r(Lya) < 0.1 where no C IV absorption is expected. 
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Fig. 8. — Median values of optical depth of C IV A1548 for different ranges of Lya 
optical depth. The open star symbols are the values measured in the HIRES spectrum 
of APM 08279+5525; no value is plotted for the bin r(Lya)= 5 — 20 because the residual 
intensity in these pixels is too low to measure the optical depth without reference to higher 
order Lyman lines. The filled circles and open squares show the results of the simulations 
described in the text. The filled circles are for the case where N(C YV)/N(H I)= 2.5 x 10~ 3 
(straight line) was assumed throughout, whereas the open squares are for simulations where 
a sudden drop by a factor of 10 in N(C IV)/iV(H I) was introduced at log iV(H I)< 14.0. 
The top and bottom panels correspond respectively to the 'ideal' and 'realistic' simulations 
discussed in the text. In the lower panel the error bars (only plotted if they are larger 
than the symbols) show the dispersion in the results when the same simulation was repeated 
twelve times. 
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Table 1. Journal of HIRES Observations of APM 08279+525 



Date 


Integration time (s) 


Wavelength range (A) 


Typical S/N 


April 1998 


1800 


4400 - 5945 


15 


April 1998 


1800 


4400 - 5945 


15 


April 1998 


1800 


4400 - 5945 


15 


April 1998 


1800 


4400 - 5945 


15 


May 1998 


2700 


4410 - 5950 


25 


May 1998 


2700 


4400 - 5945 


25 


May 1998 


900 


5440 - 7900 


30 


May 1998 


3000 


5440 - 7900 


60 


May 1998 


3000 


5440 - 7900 


60 


May 1998 


3000 


5475 - 7830 


55 


May 1998 


3000 


5475 - 7830 


55 


May 1998 


3000 


6765 - 9150 


50 


May 1998 


3000 


6850 - 9250 


50 


Summed Total 


31500 


4400 - 9250 


30 -150 
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Table 2. Details of Absorption Components Fitted to Each C IV 

System. 



System No. Redshift log N(C IV) b (km s" 1 ) 



Comments 



CI 


3.10742 


12.50 


17.2 




3.10769 


12.89 


6.7 




3.10791 


12.62 


10.5 




3.10841 


12.50 


3.7 




3.10850 


13.74 


25.7 




3.10952 


12.67 


9.2 




3.10966 


13.21 


28.2 




3.11027 


12.78 


11.8 


C2 


3.12252 


12.50 


7.8 


C3 


3.13296 


12.98 


39.0 




3.13363 


12.30 


22.5 




3.13367 


12.97 


12.4 




3.13416 


13.23 


14.9 




3.13505 


12.33 


9.6 


C4 


3.16966 


13.35 


38.5 




3.16976 


12.90 


16.3 




3.17039 


12.69 


9.8 




3.17117 


13.24 


58.6 


C5 


3.17205 


12.13 


5.7 




3.17233 


13.18 


6.7 




3.17237 


12.34 


15.8 




3.17264 


12.34 


15.8 




3.17330 


12.19 


13.3 


C6 


3.20266 


12.23 


9.9 




3.20300 


12.76 


8.8 




3.20345 


13.55 


58.9 




3.20359 


12.84 


15.8 


C7 


3.23618 


12.04 


8.4 


C8 


3.23895 


12.98 


19.4 




3.23950 


12.58 


10.5 


C9 


3.37677 


12.24 


24.5 




3.37757 


12.81 


6.5 




3.37770 


13.28 


19.7 




3.37882 


12.37 


10.8 




3.37891 


13.24 


22.5 




3.37969 


13.56 


45.4 




3.37969 


12.89 


17.0 


CIO 


3.38456 


13.13 


36.0 




3.38543 


12.65 


7.4 




3.38616 


13.51 


25.9 




3.38660 


13.02 


13.4 



log N(H I)< 14.5 



Some of 1548 A lines blended 



Some blending in both doublet lines 



1548 A line is blended with the 1550 A 
line from the previous system. 



-31 - 



Table 2 — Continued 



OjoLclIl 1NU. 




lrw N(C TVl 


h (km s~ ls ! 




Cll 


3.41088 


11.90 


1.5 


1550 A line is blended. 




3.41111 


12.35 


11.9 




C12 


3.50125 


13.01 


16.8 


Some blending in both doublet lines 




3.50141 


12.39 


9.9 






3.50204 


12.24 


5.3 






3.50209 


13.08 


22.6 






3.50281 


12.39 


32.8 




C13 


3.51380 


12.88 


17.1 






3.51436 


12.54 


14.8 




C14 


3.52706 


12.24 


13.1 






3.52760 


12.33 


16.9 




C15 


3.53303 


12.37 


8.9 




C16 


3.53744 


11.63 


2.4 


Very weak, 1550 A line barely visible 


C17 


3.55811 


12.85 


11.5 






3.55842 


12.63 


12.7 




C18 


3.57168 


11.79 


6.6 




C19 


3.58630 


12.02 


5.5 




C20 


3.62314 


11.97 


9.1 






3.62339 


11.76 


5.4 




C21 


3.63954 


12.47 


3.0 


1550 A line is blended 




3.63968 


12.09 


1.9 




C22 


3.65458 


12.90 


14.3 


1548 A line is blended 




3.65514 


12.89 


13.9 




C23 


3.66863 


12.88 


89.3 


Satisfactory fit, although slight blending 




3.66892 


13.22 


36.9 






3.67082 


12.97 


21.4 






3.67131 


12.81 


15.3 
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Table 3. Parameter Fits for Absorption Lines Shown in Figure 2 



Line No. 


Redshift 


log N(R I) 


b (km s x ) 


1 


3.34577 


14.18 


23.6 


2 


3.49289 


14.13 


25.5 


3 


3.49618 


14.28 


27.8 


4 


3.57899 


14.21 


28.3 


5 


3.58146 


14.08 


22.7 



